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Abstract 

In the case of a rarefaction fan in a non-stationary Hammersley process, we explicitly cal- 
culate the asymptotic behavior of the process as we move out along a ray, and the asymptotic 
distribution of the angle within the rarefaction fan of a second class particle and a dual second 
class particle. Furthermore, we consider a stationary Hammersley process and use the previous 
results to show that trajectories of a second class particle and a dual second class particles touch 
with probability one, and we give some information on the area enclosed by the two trajectories, 
up until the first intersection point. This is linked to the area of influence of an added Poisson 
point in the plane. 

1 Introduction 



In Hammersley (1972)} a discrete interacting particle process is introduced to study the behavior of 



the length of longest increasing subsequences of random permutations. In Aldous and Diaconis (1995), 



this discrete process is generalized to a continuous time interacting particle process on the real line, 
and they use the ergodic decomposition theorem to show local convergence to a Poisson process, 
when moving out along a ray. In this paper, we will consider Hammersley's process with sources 
and sinks, as introduced in |Groeneboom (2002) For an extensive description of this process, we 



refer to Cator and Groeneboom (2005)] since our results will be partly based on results derived 



in that paper. Here we will suffice with a brief description, based on Figure ^ 

We consider the space-time paths of particles that started on the x-axis as sources, distributed 
according to a Poisson distribution and we consider the t time axis. In the positive 

quadrant we have a Poisson process of what we call a-points (denoted in Figure H by x). At the 
time an a-point appears, the particle immediately to the right of it jumps to the location of the 
a-point. Finally, we have a Poisson process of sinks on the t-axis. Each sink makes the leftmost 
particle disappear. All three Poisson processes are assumed to be independent. To know the particle 
configuration at time s, we intersect a line at time s with the space-time paths. 

In Cator and Groeneboom (2005)1 a connection was made between the continuous time Ham- 



mersley process and the behavior of second class particles, which are well studied in the literature 



on discrete interacting particle systems such as TASEP; see for example Liggett (1999) For the 
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Figure 1: Space-time paths of the Hammersley's process, with sources and sinks. 



Hammersley process, it is natural to consider two types of second class particles: the usual one, 
where one adds an extra particle at the origin, and a dual second class particle, which corresponds 
to adding an extra sink at the origin (or removing the leftmost particle). In fact, the trajectories of 
these two particles correspond to the two longest paths of the time-reversed process such that all 
possible longest paths fall between these two longest paths. We can study the trajectories of these 
two particles at the same time, that is, for one realization of the Hammersley process. In this paper 
we study the behavior of a second class particle and its dual particle in the case of a rarefaction fan, 



a phenomenon often observed in interacting particle systems. In Ferrari and Kipnis (1995) this 



problem is considered for TASEP. In Seppalainen (2002) the Hammersley process is considered 
with general initial conditions, but the rarefaction fan is not treated. Also, our methods are quite 
different and build more on the ideas of Cator and Groeneboom (2005) In the final section we 



will study the interaction between a second class particle and its dual in the case of a stationary 
Hammersley process, and show that they will touch with probability one. We also study the area 
between the two trajectories up until this point of touch. We did not find results in the literature 
on discrete interacting particle systems that were similar to the results of our last section, so this 
interaction phenomenon may be a specific feature of the Hammersley process. 



2 Second class particles in a rarefaction fan 

Let A, fj, be two positive reals, such that \[i < 1. Let t t— > L\^(-,t) be Hammersley's process 
developing in time t, generated by a Poisson process of sources on the positive x-axis of intensity A, 
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a Poisson process of sinks on the time axis of intensity fj, and a Poisson process on of intensity 
1, where these Poisson processes are independent. Here, L\ jfM (-,t) signifies the counting process 
that counts the number of Hammersley particles on the half-line (0, do) x {t}. As was shown in 
Groeneboom (2002) the case Xfi = 1 corresponds to a stationary Hammersley process, which 



means that for each t > 0, L\„(-,t) is a Poisson process with intensity A. 

As we mentioned in the introduction, we will consider two kinds of second class particles. A 
"normal" second class particle is created by putting an extra source in the origin. A dual second 
class particle is created by putting an extra sink in the origin. The trajectories (Xt,t) of a second 
class particle and (X' t , t) of a dual second class particle are shown in Figure [2j Note that we always 
have X' t > X t . 
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Figure 2: Trajectories of (X t ,t) and {X' t ,t) 



Now consider the reversed process, where we use the North exits through [0, x] x {t} as sources, 
the East exits through {x} x [0, t] as sinks and the /3-points (these are the upper-right corners of the 
space-time paths) as our Poisson process in [0, x] x [0, t]. Burke's Theorem for Hammersley's process 



(see Cator and Groeneboom (2005) ) shows that this process is again a stationary Hammersley 
process, if we start with a stationary process. It is not hard to see from Figure|2]that the trajectories 
of X and X' correspond to longest paths from (x, t) to (0, 0) in the reversed process. 

A second class particle and a dual second class particle are symmetrical with respect to the 
main diagonal in the following sense: the reflected trajectory of a dual second class particle is a 
trajectory of a "normal" second class particle for the reflected realization, that is, with the sources 
and sinks intensity exchanged. 

A very important property of second class and dual second class particles is given in the following 
lemma, which is slightly more general than Lemma 2.2 in Cator and Groeneboom (2005)} but the 



proofs of these two lemma's are very similar, and therefore omitted here. 

Lemma 2.1 Consider a Hammersley process L with some configuration of sources and sinks. De- 
fine two coupled Hammersley processes L and L such that all three processes use the same a-points. 
Furthermore, L uses the same sources as L, but only a subset of the sinks of L, whereas L uses 
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the same sinks as L, but only a subset of the sources of L. Denote with X and X' a second class 
particle respectively a dual second class particle for the process L. Then the space-time paths of L 
and L coincide below the trajectory of X, whereas the space-time paths of L and L coincide above 
the trajectory of X' . 

The reader can verify this lemma by looking at Figure |2J 

The results in this section are also valid when A and/or p are equal to 0, if we keep in mind 
that whenever A = 0, a dual second class particle remains on the x-axis (so X' t = +00 for t > 0), 
whereas when p = 0, a second class particle remains on the t-axis (so X t = for t > 0). 

For p > 0, we denote by 1 1— > L p {-,t) a stationary Hammersley process, generated by a process 
of sources of intensity p, a process of sinks of intensity 1/p and a process of a-points of intensity 1. 

Throughout this section we are interested in the asymptotic behavior of the non-stationary 
process L\ p and trajectories of its second class and dual second class particles. We recall here the 
corresponding results for a stationary process, obtained in |Cator and Groeneboom (2005) 



Theorem 2.1 (Cator and Groeneboom) Let t — > L p (-,t) be a stationary Hammersley process 
with intensity of sources p and intensity of sinks 1/p. Let Xt be the x-coordinate of a second class 
particle at time t and X[ be the x-coordinate of a dual second class particle at time t. 
Then: 

(i) t~ 1 Xt — > 1/p 2 almost surely as t — > 00 

(ii) t X' t — ► 1/p 2 almost surely as t — > 00 



This theorem combines the results of Theorem 2.1 and Remark 2.1 of Cator and Groeneboom (2005) 
We can now show the following theorem, describing the asymptotic local intensities of the process 
L\^ when moving out along a ray t = ax, for a > 0. 

Theorem 2.2 Let a > 0. Consider the random particle configuration with counting process 

y >-> L \A X + y> ax ) ~ L *A x i ax )>y> ~ x - 

Then: 

(i) If a > 1/p 2 , the process converges in distribution, as x — ► oo, to a homogeneous Poisson process 
on R, with intensity 1/p. 

(ii) If 1/p 2 > a > A 2 , the process converges in distribution, as x — » oo, to a homogeneous Poisson 
process on R, with intensity yfa. 

(Hi) If X 2 > a, the process converges in distribution, as x — > 00, to a homogeneous Poisson process 
on R, with intensity A. 

Proof: 

(i) We couple the process L\ tfl with a stationary process -^1/^ via the same realization of the Poisson 
process in R 2 ^ , the same realization of sinks on the t-axis, and a "thickened" set of sources. This 
means that we independently add to the Poisson process of intensity A on the x-axis a Poisson 
process of intensity 1/p — A. This way of coupling (when only one set is changed (thinned or 
thickened) and the two others are kept fixed) will be used throughout the paper without further 
explanation. 

Let {X' t ,t) be the trajectory of an isolated dual second class particle at zero of the process 
As the process L\^ can be obtained from L^j^ by removing the added sources, the space-time 
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paths of both processes coincide above the trajectory (X' t ,t) (Lemma 12. Now it is enough to 
apply Theorem 12, 11 stating that X[/t almost surely converges to fi 2 , and to use the stationarity of 

£l/ M - 

(ii) This proof is analogous to the proof of Corollary 2.1, (i), i n |Cator and Groeneboom (2005) 
Let a' be such that l//i 2 > a' > a. Note that \fa' > A and l/Vff > \i. Now we couple the process 

^ with a process L generated by the same realization of points on the plane and thickened 
sets of sources and sinks of intensities \fa! and 1 / VaJ correspondingly. 

Below the trajectory (X t ,t) of an isolated second class particle at zero of L r^, the space-time 
paths of L coincide with the space-time paths of the process L which is obtained by removing 
all sinks (again Lemma l2.1|) . 

Note that for any t, the set of intersections of space-time paths of the process with the 
horizontal level t is contained in the set of intersection of space-times paths of L with the 
horizontal level t. In other words, at all times the process L ^ contains the particles of L\ tfM , and 
may have more particles, since the first process has more sources and less sinks. It follows that the 
set of intersections of space-time paths of L\„ with any horizontal interval lying to the right of the 
trajectory (X t ,t) is contained in the set of intersections of paths of L with this interval. 

Similarly, we consider a" such that a > a" > A 2 and couple with L ^tt- Now we have that 
above the trajectory (Xj., t) of an isolated dual second class particle at zero of L ^77, the space-time 
paths of L ^tjt coincide with the paths of Ll^tt, which is obtained by removing all sources. This time 
we conclude that the set of intersections of space-time paths of L\ >fl with any horizontal interval 
lying to the left of the trajectory (X' t , t) contains the set of intersections of paths of L with this 
interval. 

The situation is now as follows: if we fix R > and consider the process L\n on the horizontal 
interval [x—R, x+R] x {ax}, we know that for any fixed a' and a" such that l//i 2 > a' > a > a" > A 2 
and for x big enough, the particles of this process are a subset of the stationary process L ^ and 
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a superset of the stat ionary process L Now we ar e exactly in the situation of the proof of 
the Corollary 2.1, (i), |Cator and Groeneboom (2005~ and statement (ii) follows by considering 



the Laplace transform of the number of particles in fixed disjoint intervals and using dominated 
convergence. 

(iii) The proof of the last statement proceeds along the same lines as the proof of statement (i). 
Couple L\a with L\, note that their space-time paths coincide below the trajectory (Xt,t) of an 
isolated second class particle at zero of the process L\ and use Theorem 12 . 1 1 and the stationarity 
of L\. 

□ 

Let F\^(x,t) be the flux of the process L\ ilM at (x,t), which we define by the number of space- 
time paths of contained in the rectangle [0,x] x [0,i\. Note that F\ itl (x,t) can be represented 
as the number of crossings of the interval [0, x] x {t} by space-time paths of L\n plus the number 
of crossings of the interval {0} x [0, t]. It is also equal to the number of crossings of the interval 
[0, x] x {0} plus the number of crossings of the interval {x} x [0, t] by space-time paths of L\^. 

Definition 2.3 When moving out along a ray t = ax, we define the local intensity function I(x,h) 
by the following formula: 

I(x, h) = E(F Xjfi (x + h, ax) - F Xjfi (x, ax)). 
We suppress the dependence on a, A and [i in our notation. 
Note that 

I(x, h) = E(L\^(x + h, ax) - L x ^(x, ax)), 

so it is equal to the expected number of crossings of the interval [x, x + h] X {ax} by space-time 
paths of the process L\ t/Ji . 

Corollary 2.4 Fix the positive constants a, A and jj,. We have 

{hjii, a > 1/V 2 , 

Ja ■ h, 1/u 2 > a > A 2 

A • h, A 2 > a > 0. 

Proof: Coupling in a standard way the process with the process -^i/^ (adding sources), we 
can estimate 

L\^{x + h, ax) - L x ^{x, ax) < L l/lx (x + h, ax) - L 1/fJ _(x, ax) 
for any realization and any x > 0. Define for x > 

A x = L\^(x + h, ax) - Lx ifl (x, ax) and B x = L 1/fl (x + h, ax) - Ly^x, ax). 

Define a random variable A, whose distribution depends on the fixed parameters a, A and fi, as 
follows: 

'Pois(/i///), a > 1/fJ?, 
A ~ < Pois(^/a • h), l/ / u 2 >a>A 2 
Pois(A-/»), A 2 >a>0. 



G 



Theorem 12.21 shows that 

Furthermore, 

For any M > we have 



A x —> A (x — ► oo). 



V x > : ~ Pois(/t///). 
E(A,) = A M) + E((^ - M)l {Ax>M} ). 



Also, 

V x > : E((4« - M)l {Ax>M} ) < E((B a - M)l {Bx>M} ) = E((S - M)l {Bo>M} ). 



Since 



and 



lim E(A„ A M) = E(y4 A M), lim E(A A M) = ~EA 

z— >oo M— >oo 



lim E((B - M)l {So>M} ) = 0, 



we conclude that 

lim E^) = EA. 



□ 



In the next theorem we will discuss the behavior of Xt and X[, the second class particle and 
the dual second class particle of the process Lx^. 

Theorem 2.5 Let be the non- stationary Hammersley process defined above. Let Xt be the 
x-coordinate at time t of a second class particle starting at zero. Then 



lim P(X ax > x) 



fl, a>l/v 2 , 

tH, i/m 2 >«>a 2 

0, A 2 > a > 0. 



Proof: To prove this theorem we will calculate the derivative J^I(x, h)\h=o- The proof is inspired 
by the proof of Theorem 1, Ferrari and Kipnis (1995)[ where the analogous problem for TASEP 



is solved. 

Fix t and x and consider some small h > 0. We couple the process L\^ with a process L\ t a 
constructed in the following way. Consider the same Poisson realizations of points in the plane. 
Consider the same realization of sinks. Consider the same realization of sources on the interval 
(h,oo) x {0}. On the interval [0, h] x {0} we add an independent Poisson process of sources of 
intensity 1/fj, — A. The situation is illustrated in Figure |IJ 

Note that the process L\^, restricted to the rectangle [0, h] x [0, t], is a stationary Hammersley 
process with intensity of sources 1 j \x and intensity of sinks [i. Thus, the set of crossings of space- 
time trajectories of L\^ with the vertical segment {h} x [0,t] is a Poisson process of intensity 
H, independent of the sources on (h, oo) x {0} and of the Poisson points in (h, oo) x [0, t] (this 
follows from the stationarity of the Hammersley process when viewed from left to right). Now it 
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Figure 4: trajectory Xj 1 of added point (•) hits the right side and the flux F\ )tl (x,t) is decreased; 
for this realization it decreases from 4 to 3. 



follows that the number of space-time paths of the process L\ tfJi inside the rectangle [h, x + h] x [0, t] 
(which we will denote by F\^(x,t)), has the same distribution as the flux of the process on 
the rectangle [0,x] x [0,t]. 

We are interested in the difference 

F x ^(x + h,t)-F Xifl (x,t). (2.1) 

This difference is equal to the difference between the number of space-time paths of the process 
in the rectangle [0, x + h] x [0, t] and the number of space-time paths of the process in 
the rectangle [h, x + h] X [0, t\. 

As we are interested in the expectation of this difference only up to the first order of h when 
h I 0, we can distinguish two cases. In the first case there is no extra source in [0, h] x {0}, which 
means that the space-time paths of ^ and ^ are the same. In the second case there is an 
extra source in [0, h] x {0}. 

In the first case the difference (|2.1(l is equal to the number of sources of in [0, h] x {0}, since 
those sources are not counted in F\^(x,t). Clearly, the expected number of sources in [0, h] x {0} 
equals Xh. 

In the second case we have to be more careful. First of all, the probability of having an extra 
source equals, in first order, (l//i — X)h. The added source can affect the space-time paths of the 
original L\ „ inside the rectangle [h, x + h]x [0, t] and can change the flux. Note that the change in 
the space-time paths will happen above the trajectory of the added point, considered as a second- 
class particle. If this trajectory hits the right side of the rectangle (i.e. {x + h} x [0, t]), it follows 
that one of the intersections of space-time paths of with the right side is eliminated in L(A, /j,); 
see Figure |IJ Thus, in this case the difference (j2.1|) is +1. But if the trajectory of the added point 
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(which we define by X^ ) does not hit the vertical segment {x + h} x [0,t], it follows that the sets 
of intersections are the same for both processes and the difference is zero again. 
Now we can conclude that 

E(F X!tl (x + h,t)- F x ^(x, t)) = \h + - \)hP(X? >x + h)+ 0(h 2 ). (2.2) 

Taking t = ax, dividing by h and taking the limit h — > we have 

J-I(x, h)\ h=0 = A + (l//i - A)P(X a;E > x). 

Prom this it follows that is non- increasing in h € [— 1, 1] for any fixed x > 1. Indeed, fix 

ho G [-1,1]- Use (f2~2~)) to conclude that 

d d 

— I(x, h)\ h =h = q^I(x + ho, h)\ h=0 = A + - A)P(X a:r > x + h Q ). 

The probability P(X ax > x + ho) obviously does not increase in ho. 

We are now in the following situation: we have a family {I(x, h) : x > 1} of functions of h, where 
h £ [—1,1]. For all x the function I(x,h) is differentiable on [—1,1] and concave. Moreover, we 
know that the family I(x, h) has a pointwise limit, which is also differentiable on [—1, 1] (Corollary 
12. 4|) . It follows that the family of derivatives also converges to the derivative of the limit, in other 
words 

1 1 

lim lim —I(x,h) = lim lim —I(x,h). 

x^oo h-+0 h h— >Q x^oo h 

This is a standard result from analysis. 
It follows that 

1/fX, a > 1/V 2 , 

y/a, 1/fi 2 > a > A 2 

A, A 2 > a > 

which finishes the proof of the theorem. 

□ 



A + (l//i - A) lim P(X ax >x) = { 

x—>oo 



Corollary 2.6 Let L\„ be the non- stationary Hammersley process defined above. Let X[ be the 
x-coordinate at time t of a dual isolated second class particle at zero. Then 



lim P{X' ax > x) = < 



1, _ 

x/x-JTfc 



l/A-At, 



0. 



a > l/n 2 , 
1/H 2 > a > A 2 
A 2 > a > 0. 



Proof: This follows from the symmetry of the trajectories of a second class particle and a dual 
second class particle (reversing the role of x and t, interchanging A and [i and replacing a by 1/a). 

□ 
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3 Second class particles and dual second class particles in the 
stationary case 

Now we will turn our attention to the stationary case. Let A > and consider the Hammersley 
process L\ with intensity of sources A, intensity of sinks 1/A and intensity of a-points 1. As it 
follows from Cator and Groeneboom (2005)1 ( see a ^ s0 Theorem 12.1(1 . the trajectories of a second 



class particle at zero, (X t ,t), and of a dual second class particle, (X' t ,t), both asymptotically tend 
to the ray t = X 2 x. We will prove that these trajectories intersect with probability one and provide 
some information on the area enclosed by these trajectories, the origin and the first intersection 
point. 

Theorem 3.1 Let L\ be the stationary Hammersley process defined above and let (Xt, t) and (X' t , t) 
be the trajectories of a second class and a dual second class particle respectively. Then (X t ,t) and 
(X' t ,t) intersect with probability one. 

Proof: To prove this theorem we will again use coupling. Consider the standard coupling of L\ 
with L\fi by removing all sinks. The space-time paths of L\ and La,o coincide strictly below the 
trajectory (Xt,t) (Lemma 12.1(1 . Choose a realization such that (X^,t) tends to the ray t = X 2 x. 
If (X[,t) and (Xt,t) do not intersect for this realization, then the trajectory of a dual second 
class particle of the process L\ q coincides with (X' t ,t) and, consequently, will also tend to the ray 
t = X 2 x. Suppose this event has probability e > 0. Denote X^t) as the x-coordinate at time t of 
a dual second class particle for the process L\ . The previous discussion shows that 

Pi lim X' (X 2 x)/x = !)>£. 

x— >oo 

Then for any a > A 2 , we would have 

A 2 

P(]M:Vi>M: X' G {X 2 x) > —x) > e 
This implies that there exists M > such that 

A 2 



P(V x > M : X' (X 2 x) > —a?) > e/2. 



This in turn implies 



A 2 a 
V x > M : P(X' Q (X 2 x) > — x) > e/2 O V y > -=M : P(X' (ay) > y) > e/2 

A 



This contradicts Corollary 12.61 for a close enough to A 2 , which proves that (Xt,t) and (X' t ,t) inter- 
sect with probability one. 

□ 



Now we will study the expected area enclosed by the trajectories of a second class and a dual 
second class particle. The main idea comes from the following observation: suppose we have a 
realization of a stationary Hammersley process in the positive quadrant. Now we add an extra 
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Figure 5: In the left picture, the flux at the point (x, t) increases, for this realization from 2 to 3. 
In the right picture, the flux at (x, t) remains the same, for this realization equal to 2. 



a-point at (u, v). How does this change the space-time paths? The situation is depicted in Figure 

m 

Adding an extra a-point at (u, v) moves the first particle immediately to the right of (u, v) to 
(u,v); this corresponds to a dual second class particle that starts at (u,v). However, since there 
is now a particle at (u,v), the first "sink", or crossing on {u} x (v,oo), also disappears, which 
corresponds to a "normal" second class particle that starts at (u,v). Within the region, enclosed 
by the trajectories of the second class particle and the dual second class particle, the space-time 
paths will change in such a way that the flux F(x, t) is increased by 1 for each (x, t) within the 
trajectories. However, the first time the two trajectories meet (which happens with probability 1, 
according to Theorem I3.1|) . the two particles annihilate each other (since one corresponds to an 
extra source, whereas the other corresponds to an extra sink). This means that the shaded area in 
Figure [SJ between the two trajectories, is the only area where the space-time paths are changed, 
and within this area the flux is increased by 1. 

Now we define the following "rainfall" process, which will emulate the adding of extra a-points 
to a stationary Hammersley process. We define this process in [0, oo) 3 , where the first two variables 
x and t are the same as we have seen before, but the third parameter r is a kind of external 
time parameter. Consider three independent Poisson processes in this space: a Poisson process 
of intensity 1 on the plane [0, oo) x {0} x [0,oo), a Poisson process of intensity 1 on the plane 
{0} x [0, oo) x [0, oo) and a Poisson process with intensity measure 2rdxdtdr on [0, oo) 3 . 

Then define as the Hammersley process at "time" r that uses as its sources the projection 
of all Poisson points in [0, oo) x {0} x [0, r] onto the x-axis, as its sinks the projection of all Poisson 
points in {0} x [0, oo) x [0, r] onto the t-axis and as its a-points, the projection of all Poisson 
points in [0, oo) 2 x [0, r] onto the (x, t)-plane. It is easy to check that will be a stationary 
Hammersley process with source and sink intensity equal to r, and the intensity of a-points equal 
to r 2 . The advantage of this construction lies in the fact that we have coupled a family of stationary 
Hammersley processes. 

We will use this process to prove the following theorem. 

Theorem 3.2 Consider a Hammersley process with intensity of sources and sinks being r and 
intensity of a-points being t 2 . Let be the region in the plane enclosed by the trajectories 
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(X T (t),t) and (X' T (t),t) of a second class and a dual second class particle starting in the origin, 
until their first intersection point. Fix x,t > and let Area denote Lebesgue's measure on the 
plane. Then 

EArea(0 (r) n [0,x] x [0,t]) = - — (E(x - X T (t))+ + E(t - X T (x))+). 

Proof: Consider the rainfall process described above and let F T (x,t) be the flux at the point (x,t) 
of the stationary process L( T ). To prove the theorem we will calculate the derivative 

^-MF T (x,t) 

OT 

in two different ways. 

Take some small dr > and note that, since the processes and L^ T+dr ^ are stationary, we 
can calculate the expected flux at both times directly: 

EF T+dr (x, t) - EF T (x, t) = x(t + dr) + t(r + dr) - xr - tr = (x + t)dr 

and therefore we get 

d 

—EF T (x,t) = x + t. (3.3) 

OT 

Now we consider a different way of calculating the difference 

E(F T+dr (x,t) -F T (x,t)). (3.4) 

by using the coupling of and L^ T+dT \ Up to the first order of dr we can assume that in the 
time interval [r, r + dr] there either falls one extra source onto the segment [0, x] x {0}, one extra 
sink onto the segment {0} x [0, t] or one extra a-point inside the rectangle (0, x] x (0, t] (if nothing 
happens the difference is of course zero). 

In the first case the flux will only be changed if the trajectory of a second class particle, starting 
at the dropped source on [0,x] x {0}, will cross the line [0, oo) x {t} to the left of the point (x,t). 
In this case the flux will be increased by one. Similarly, in the second case the flux will increase by 
one only if the trajectory of a dual second class particle starting at the dropped sink on {0} x [0, t] 
will cross the line {x} x [0, oo) to the right of the point (x,t). 

The area we are interested in arises from the third case. It follows from the explanation around 
Figure El that the flux at (x,t) is increased by one only if (x,t) £ 0^ T \u,v), where (u, v) is the 
location of the extra a-point and 0^ T \u,v) denotes the area described in the theorem, but here 
the second class particle and the dual second class particle start at (it, v). 

Combining all three cases we have 

X t 

E(F T+dT (x,t) - F T (x,t)) = dr J P(X Ttt (u,0) < x)du + dr J P(X' rt (0,v) > x)dv + 



(3.5) 

X t 

+ 2rdr j j P({x,t) £ {r) (u,v))dudv + o(dr), 
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where X Tjt (u, 0) and X' T t (0, v) are x-coordinates at time t of a second class and a dual second class 
particle of L^ T \ starting at points (u, 0) and (0, t>) respectively. 

Using the fact that in the stationary Hammersley process LW, OW and 0^ T \u,v) have the 
same distribution modulo a translation, we can rewrite the last integral of ()3.5|) in the following 
way: 

X t X t 



= E y y l oW (u,t;)dtHfo = EArea(0 (r) n [0,x] x [0,t]). 
o o 

Dividing both parts of 1)3. 5 Jl by dr, taking the limit dr — > and using 1)3.3)) we arrive at 



x 



x i 

+ t = y P(X T>t (u, 0) < x)du + y P(^ )t (0, u) > s)d« + 2rEArea(0 (r) n [0, x] x [0, t]). 



Again using the stationarity of , we can see that 

P(X T ,t(ii, 0) <x) = P(X T (t) <x-u) 

and 

P{X' Ttt (p,v)>x) = P(X' T (t-v)>x). 

This leads to 

X t 

EArea(O r n [0,x] x [0,t]) = - — / P(X T (t) < u)du - — f P(X' T (v) > x)dv. (3.6) 

2r 2t J 2t J 

o o 

Consider the first integral on the right side of 1)3.6)1 . We can rewrite it in the following way: 

X 

P(X T (t) < u)du = nj l{x T (t)<u}du) 



= E(x-X T (t))+. 

Using the equality P{X' T {v) > x) = P(X T (x) < v), which follows from the symmetry of the process 
L( T ) (note that we have the same intensity of sources and sinks), we conclude a similar formula for 
the second integral in 1)3.6)1 



y P(X' T (v) > x)dv = E(t - X T (x)). 



o 

Combining these results we have proved that 

EArea(0 (T) n [0, x] x [0, t]) = - — (E(x - X T (t))+ + E(t - A T (x))+). □ 

2r 2r 
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This result is only valid for stationary Hammersley processes with equal source and sink in- 
tensity. However, if we consider L x , a stationary Hammersley process with source intensity A, 
sink intensity 1/A and intensity of a-points equal to 1, it is clear that this process has the same 
distribution as the image of the process L\=i under the map 

$\(x,t) = {x/X,Xt). 

This observation easily leads to the following corollary. 

Corollary 3.3 Consider a stationary Hammersley L x . Let O x be the area in the plane enclosed 
by the trajectories (Xf,t) and (X[,t) of a second class and a dual second class particle starting in 
the origin, until their first intersection point. Fix x,t > and let Area denote Lebesgue's measure 
on the plane. Then 

EArea(0 A n [0, x] x [0, *]) = I 1 - -(E(x - X t )+ + E(i/A 2 - X x , x ) + ). 

Another interesting asymptotic corollary to Theorem 13.21 is the following. 
Corollary 3.4 Consider a stationary Hammersley process L x . 

1. For any fixed xq > we have 

Axn + AEX\2„„ 
EArea(0 A n [0, x ] x [0, oo)) = — 

2. For any fixed to > we have 

EArea(0 A n [0,oo) x [0,t ]) 



t /X + XEX t0 
2 



EArea(0 A n [0, x] x [0, X 2 x\) _ 
lim — A. 

x^oo x 

Proof: For the first statement, consider E(x — Xt)+. Note that 

(x — Xt)+ — » a.s. 

as t — > oo. Since (x — X t )+ is bounded by x, we get 

E(x-X t ) + ^0 (i -►(»). 

Now consider t/X 2 - E(t/X 2 - X X 2 X ) + = E(mm(X X 2 x ,t/ X 2 )). Clearly, 

E(min(X A2;E ,t/A 2 )) -» EX X 2 X (t -» oo), 

so the first statement is a direct consequence of these observations and Corollary 13.31 The proof of 
the second statement is completely analogous. 

The last statement follows from the fact that in L x , we have 

fe^ - 1 a.s. 
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Now we use dominated convergence to get the desired result from Corollary 13.31 



□ 



This last corollary seems to indicate that the two trajectories of the second class particle and 
the dual second class particle run parallel at average vertical distance A. However, the true behavior 
of the two trajectories is that they often collide quite quickly, but with low probability they drift 
away from each other, in which case 0\ n [0, x] x [0, X 2 x] grows quadratically in x. It does follow 
from Theorem 13. 1 1 that eventually the two trajectories will touch. We wish to remark that it is not 
hard, using the methods of Cator and Groeneboom (2005) and this paper, together with some 
crude estimates, to show that EX t < +oo. 
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